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The aim of this article is to present a time-frequency theory for orthogonal 
polynomials on the interval [—1,1] that runs parallel to the time- frequency 
analysis of bandlimited functions developed by Landau, Pollak and Slepian. 
For this purpose, the spectral decomposition of a particular compact time- 

i 1 frequency-operator is studied. This decomposition and its eigenvalues are 

closely related to the theory of orthogonal polynomials. Results from both 
theories, the theory of orthogonal polynomials and the Landau-Pollak-Slepian 
<-h theory, can be used to prove localization and approximation properties of the 

corresponding eigenfunctions. Finally, an uncertainty principle is proven that 
reflects the limitation of coupled time and frequency locatability. 
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2 1 Introduction 

In the beginning of the 1960s, Landau, Pollak and Slepian developed a remarkable theory 
on the time- frequency analysis of band-limited functions. In a series of papers (|22|. |23j . 
|24j . |37| . |38j . |40| ) they studied the interplay between the two projection operators Pa 
and Pg defined on the Hilbert space L 2 (R) for two intervals 4,Bclby 

P A f:=XAf, PTf:=XBf, feL 2 (R). 
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1 Introduction 



They analyzed the composition PbPaPb an d its spectrum and found that the eigenfunc- 
tions of the compact self-adjoint operator PbPaPb are well-known special functions: the 
prolate spheroidal wave functions. Using these particular eigenfunctions as a basis for 
the band-limited functions in L 2 (M) on the other hand, they were able to prove a series 
of interesting results concerning the approximate concentration of functions in the time 
and the frequency domain, as well as an uncertainty principle involving a lower bound 
for the angle between the vectors P^f and PbJ ■ An overview of these results can be 
found in the articles |25j . [39j and the book [5, Section 2.9]. 

Later on, the Landau-Pollak-Slepian-theory was extended to a variety of different set- 
tings. Among others, there exist analogies on the unit circle [39], on discrete groups 
|18] and on symmetric spaces like the unit sphere |19j . |36j . Various generalizations of 
this theory can be formulated, for instance by considering eigenfunctions of particular 
Sturm-Liouville differential equations [42j or using reproducing kernel Hilbert spaces |43| . 
Particularly interesting for this article is the fact that there exists also an extension of 
this theory to orthogonal polynomials defined on subsets of the real line [32J. 



The aim of this paper is to present a time-frequency analysis for orthogonal polynomials 
on the interval [—1,1] that runs parallel to the Landau- Pollak-Slepian theory described in 
|32j . For the frequency localization of a function / in the weighted L 2 -space L 2 ([— 1, 1], w) 
we will use, as in [32], an operator P™ that projects the function / onto a finite dimen- 
sional polynomial space II™. However, in contrast to the theory outlined above, we will 
not use a projection operator Pa to describe the space localization of /. Instead, we will 
consider the multiplication operator M x defined by multiplying the function / with the 
variable x. 

Compared to the projection operator Pa, the usage of the multiplication operator M x 
leads to a time-frequency analysis in which the localization of / at the boundary points 
x = 1 and x = — 1 of the interval [—1,1] plays an important role. For a normalized 
function / G L 2 ([— 1, 1], w), the mean value £(/) — (M x f,f) w is located in the interval 
(— 1, 1). The closer e(/) gets to 1 or —1, the more the L 2 -mass of / is concentrated at 
x = 1 or x = — 1, respectively. Therefore, the mean value e(f) can be considered as a 
measure on how well the function / is localized at the boundary points x = 1 or x = — 1. 
Particularly this property of e(/) implies the possibility to construct polynomials in II m 
that are optimally localized at the boundary of [—1, 1] (see [7], |16j . [33]). 

The principal examination object for the time-frequency analysis in this paper is the 
finite dimensional self-adjoint operator P™M X P™ in combination with its eigenvalues 
£ m fc) 1 < k < n — m + 1, and corresponding eigenfunctions ip™ k - One of the main 
advantages of the operator M x in place of Pa is the fact that the spectral decomposition 
of P™M X P™ is closely linked to the theory of orthogonal polynomials. This relation 
makes it possible to use a very large repertoire of techniques and results from the theory 
of orthogonal polynomials to analyse the properties of the spectral decomposition of 



P™M X P™. In the spectral Theorem 2.1 we will see that the eigenvalues of P^M X P^ 
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2 The spectral decomposition of P™M X P™ 
are precisely the roots of the associated orthogonal polynomials Pn—m+i 

(x,m). Also 

the eigenfunctions can be stated explicitly. In the case m = 0, they correspond to the 
fundamental polynomials of Lagrange interpolation. 

A second advantage of using the operator M x consists in the fact that the value e(f) 
represents also the expectation value of the L 2 -density /. The density / can be considered 
as localized at the expected value e(f) if the variance var(/) is small. Therefore, we 
can investigate the localization properties of the eigenfunctions t/j™k °^ F^-M^P™ by 
considering the variances var(^, m fc ). In order to show that the functionals var(^™ fc ) are 
small when n is large, we will use results of Nevai, Zhang and Totik [29] . |30| on uniform 
subexponential growth. The major result in this context is Theorem |3.6| It states that if 
the weight function w of the space L 2 ([— 1, l],w) is in a particular subclass of the Nevai 
class M(0, 1), then the variance of the eigenfunctions ip™ k tend to zero as n — > oo. 

In Section |4j we will analyse how the decomposition of a bandlimited function / G II™ in 
the single eigenfunctions il>™k can ^ e use( ^ to approximate functions that are localized at 
a point or a subinterval of [—1, 1]. In this case, not all the eigenfunctions ift™ k are needed 
to approximate the function /, but just those that are situated in the region in which 



/ is concentrated. In Theorem 4.1 and 4.3 we will give simple error estimates for such 



approximations if the function / is localized in a certain area or at a particular point of 
the interval [—1,1], respectively. 

Finally, we will prove an uncertainty relation for orthogonal polynomials involving the 
operators M x and P™. This relation can be considered as an extension of the angular 
uncertainty principle in the Landau-Pollak-Slepian theory. For a normalized function 
/ C L 2 ([— 1, l],u>), the determining quantities of the uncertainty relation are the norm 
ll-PnVllw aR d again the mean value e(/). The norm H-F^/H™ gives a measure on how 
well the function / is concentrated in the polynomial subspace II™. On the other hand 
the value e(/) can be seen as a measure of the localization of / at the boundary points 
x = — 1 and x = 1. The main result in the last section is Theorem |5.5| claiming that for a 
normalized function / G L 2 ([— 1, l],w) it is impossible that and are both 

close to 1. In particular, this result implies that if |e(/)| is too close to 1, / cannot be a 
polynomial in II™. 



2 The spectral decomposition of P™M X P™ 



We consider the Hilbert space L 2 ([— 1, l),w) with the inner product 



(f,g)w ■= j f(x)g(x)w(x)dx, 

and a positive weight function w having finite moments J_ x n w(x)dx, n G N. By {p/}^ , 
we denote the family of polynomials pi of degree I that are orthonormal on [—1, 1] with 
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respect to the inner product (•, Further, we assume that the polynomials pi are 
normalized such that the coefficient of the monomial x l is positive. Then, the family 
{pi}lZ defines a complete orthonormal set in the Hilbert space L 2 ([— 1, l],w) (cf. |41| 
Section 2.2]). By U n , we denote the polynomial space spanned by the polynomials p\ up 
to degree n, and by n™ the polynomial wavelet space spanned by the polynomials pi, 
m < I < n. 

For a normalized function / G L 2 ([— 1, l],iu), \\f\\ w = 1> we define the mean value e(f) 
and the variance var(/) by 

e(/) := ^ x\f(x)\ 2 w(x)dx, (1) 
l 

2| f/„M2.,.c \ t / _2i f/ m2„ 



var(/):=/ (x - e(/))V(x)|Mx)ds = / ^|/(x)|^(x)da: - (2) 



i 



/i X2|/(: 



We are now going to introduce a time-frequency analysis for functions / £ L 2 ([—l,l],w) 
based on the following two operators: 

(M x f)(x):=xf(x), (3) 
n 

(P™f)(x):=Y,(f,PiUPi(x)- (4) 

l=m 

If m = 0, we write P n instead of P^. The multiplication operator M x as well as the 
orthogonal projection P™ onto II™ are both self-adjoint and bounded operators on the 
Hilbert space L 2 ([— 1, 1], w). Therefore, also the composition 

P?M X P™ (5) 

is a bounded and self-adjoint operator on L 2 ([— 1, 1], w). Moreover, since P™ is compact, 
P™M X P™ is also a compact operator. Hence, by the Hilbert-Schmidt theorem the spec- 
trum of the operator P^M^P™ is discrete (it is even finite) and the eigenfunctions form 
an orthogonal basis of L 2 ([—l, 1], w) (cf. |34[ Theorem VI. 16]). The subsequent Theorem 



2.1 will illustrate that the eigenvalues and eigenfunctions of P™M X P™ are well-known 



in the literature. 



For a description of the spectral decompostion of P™M X P™, we need first of all the 
notion of associated polynomials. We know that the orthonormal polynomials pi satisfy 
the three-term recurrence relation (cf. |17[ Section 1.3.2]) 

k+lPl+l(x) = (x - ai)pi(x) - bipi-x(x), I = 0,1,2,3,... (6) 

p-i{x) = 0, p (x) = —, 

with coefficients o; S R and 6/ > 0. For m G N, the associated polynomials pi(x,m) on 
the interval [—1,1] are then defined by the shifted recurrence relation (see |17[ Section 
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1.3.4], [2D Section 2.10]) 

b m+ i +1 pi + i(x,m) = (x - a m+ i)pi(x,m) - b m +iPi-i(.x,m), 1 = 0,1,2,..., (7) 
P-i(x,m) = 0, po(x,m) = 1. 

For m = 0, we have the identity, pi(x, 0) = boPi(x). The polynomials p;(x) and pi(x, m) 
can be described with help of the symmetric Jacobi matrix J™ < m < n, defined by 



J m+1 







-1 a m+l m+ 2 










V o 



b m +2 O-m+2 b m+ 3 





o \ 








bn-2 On-1 &n-l 

6 n _i a n / 



(8) 



If m = 0, we write J n instead of J^. Then, in view of the three-term recurrence formulas 
([7]), the polynomials p\ and pi(x,m) can be written as (cf. )21[ Theorem 2.2.4]) 

, , 1 



piyx) 
Pi[x,m) 



bo 



det(xl; - J|_i), 



det(slj- «C+/_i), 



(9) 
(10) 



where 1/ denotes the ^-dimensional identity matrix. We can now explicitly state the 
spectral decomposition of the operator P™M X P™ . 

Theorem 2.1. The operator P™M X P™ on L 2 ([— 1, 1], w) has the spectral decomposition 

n—m+l 



PnM x P™f— x n,kUii : n ) k)w' t P™k- 



fc=i 



For m > 1, the eigenvalues x™ k denote the n — m + l roots of the associated polynomial 
p n _ m +i(x,m) and the eigenf unctions have the explicit form 



^n,k\ x ) — K n,k 



with the normalizing constant 



1 p n+1 (x)p 

n—m 



X — X 



v n,fc 



PJ-m«' 



n,k 



™ k ,m)^ 2 



(12) 



(13) 



l=m 



For m = 0, the eigenvalues x n \. correspond to the n + 1 roots of the polynomial p n+ i(x) 
and the eigenf unctions ip n k correspond, up to a normalizing factor, to the fundamental 
polynomials of Lagrange interpolation, i.e. 



1pn,k( X ) = K n,kPn(x n< k)bn+ 



Pn+l(x) 
X X n h 



(14) 
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where 



1 



K n,k ■= (^2pi(x n ,k) 2 ) 2 ■ (15) 



1=0 



Proof. We consider the projection P™f of the function / onto the subspace IL™ in terms 
of the expansion P™f = s ^d=m c ^'Pl with the coefficients c\ = (f,pi) w - Using the three 
term recurrence relation (J7| it is straightforward to show (see [7J Lemma 2.7]) that the 
mean value e(P™f) of P™ f can be written as 



(P™M x P™f, f) w = (M x P™f, P™f) w = e(P™f) = c" J-c, (16) 

where c H denotes the conjugate transpose of the vector c = (c m , . . . ,c n ) T . Thus, the 
eigenvalues of P™M X P™ in II™ C L 2 ([— 1, 1], w) correspond to the eigenvalues of the 



Jacobi matrix J™. On the other hand, by equation (10) the eigenvalues of J™ are 
exactly the roots x™ k , k = 1, . . . n — m + 1, of the associated polynomial p n _ m+ i(x, m). 
The eigenvector c& corresponding to the root x™ k is simple and can be computed via the 
three-term recursion formula ([7]) as 

Cfc = (l,pi(x™ k , m),...,p n - m (x™ k ,m)) T . (17) 

The corresponding normalized eigenfunction ip™ k of P™M X P™ can then be written as 

n 

i>n,ki x ) = <fc Yl Pl-™( x n,k, m) P l(x), (18) 
l=m 

with the normalizing constant given in ( |13| ). By an alteration of the classical 

Christoffel-Darboux formula (see |7J Lemma 3.1]), the eigenfunctions for m > 1 
have the explicit form 

lPn+l(x)p n—m 

Vn,k\ X > = K n,k _ m ■ 

For m = 0, we get directly by the Christoffel-Darboux formula (see [H Chapter 1, 
Theorem 4.5]) that 

Wn,k\ x ) = K n,fc0n+1 ~ • 

□ 



Remark 2.2. In the literature, the spectral Theorem 2.1 is well-known for the case m = 
(cf. |1] Lemma 8.4] and |35| Proposition 1.3.1]). For the more general case m > 0, an 
equivalent representation of Theorem 12. II is the eigenvalue decomposition J™c^ = £™ fc Cfc 



of the matrix J™ (see |17[ Section 1.3]) . lb the best of the authors knowledge, the explicit 



formulas (12) of the eigenfunctions ip™ k , m — 1> can ^ e considered as novel. 
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Remark 2.3. The eigenfunctions {V^fcifc^i"^ 1 °f the operator P™M X P™ form an or- 
thonormal basis of the polynomial space n m . Hence, we can expand polynomials P £ II™ 

as 



n— m+1 

p{x)= £ (p,w,kU<k(x)- 

k=l 



In the case m = the functions ip n k correspond to the fundamental polynomials of 
Lagrange interpolation and can be described through the Christoffel-Darboux kernel (see 



|27[ (1.1.9)] and formula (18)). The functions ip n ,k a re used in [12] and [13] as particular 
orthogonal scaling functions in a wavelet decomposition of a function / £ L 2 ([— 1, 1], w). 
If m > 1, the construction of the wavelet basis functions in these two papers differs 
however from the eigenfunctions ip™L. considered in this article. For a general overview 
on polynomial frames and polynomial wavelet decompositions, we further refer to the 
articles HQ, [28] and the book [27] . 

Remark 2.4. It was specified in the introduction that the mean value s(f) can be inter- 
preted as a measure on how localized the function / is on the boundary points x = 1 and 
x = —1 of [—1, 1]. In the following, we will say that a function / is localized at x = 1 
or x = —1 if the mean value e{f) approaches 1 or —1, respectively. For a polynomial 
P £ IT™, the mean value e{P) can be written as e(P) = (P™M X P™P, P) w . Precisely 
this mean value e(P) was used in |6] and [7] to construct polynomials in II n and Il m that 
are optimally space localized at the boundary points x = 1 and x = — 1 of the interval 
[—1,1]. These optimal polynomials are exactly the eigenfunctions V'rTmax an d ^rTmin m 



Theorem 2.1 corresponding to the largest and the smallest eigenvalue of the operator 



P™M X P™. By (16), we have for the largest eigenvalue of P™M X P™ the relation 



max (P™M X P™P, P) w = max c H J m c. 

^en™,|[P|| TO =i c H c=i 



This characterization is thoroughly used in |16j to get estimates for the largest zero of 
orthogonal polynomials. 

Taking a step further, we can also consider the orthogonal complement II m 0span{^™ max } 



of ipn'mzx in II™. Then, the spectral Theorem 2.1 says that the polynomial in II™ 
span{-(/>™ max } that is best localized at x = 1 is the eigenfunction V'rTmax-i corresponding 
to the second largest eigenvalue a; m max _i of P™M X P™ . Hence, repeating this argumen- 
tation, Theorem 



which the fc-th. 



2.1 produces a chain of elementary orthonomal basis functions in 
element is worse concentrated at x = 1 than the (k + l)-th. element 
V'JTfc+i but better than the (k — l)-th. element •i/Cfc-r ^he measure of the corresponding 
localization is given by the mean value ^V'JTfc) = x n it- 
Example 2.5. We consider the orthonormal Chebyshev polynomials t n of first kind defined 
by (see PS p. 28-29]) 

io(cost) = t n (cost) = Y^fcos(nt), n > 1, cost = x. 

The roots of the Chebyshev polynomials t n +i are given by x n ^ = cos 2 " 2 ^ 2 | fc ^ l ~ 3 7r, k = 
1, . . . ,n + 1 (see |41[ (6.3.5)]). The normalized associated polynomials t n (x,m), rn > 1. 
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correspond to the Chebyshev polynomials u n of the second kind given by (see |17[ p. 28- 
29]) 

sin(n + l)t 



-u n (cos t) 



sini 



n > 0. 



The zeros of the polynomial u n _ m+ i are given by x™ k = cos n ~ r ^^ L ^J 2 k 7r , k = 1, . . . , n — 
m + 1. Hence, by the formulas (12) and (14) we get for the eigenfunctions -0™ fc the 
following explicit representation 



V>n,fc(cOsf) 



Kn,k COS 2n+2 



n(2n-2fc+3),r ^ + ^ 



vr cos t - cos 2n ^ 2k + 3 7r 



2n+2 



-1) 



n— m— k+1 



:+ cos(n + l)t + cos(m — l)i 



2vr 



cost- cos ™^vr 



m > 1. 



The constants K n & can be computed explicitly and are given as (see [27, Formula 
(1.1.17)]) 

/ _ 2ra + 1 + U 2 n{x n ,k) 

Some of the eigenfunctions tl>™k are illustrated in Figure [IJ 
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In this section, we are going to investigate localization properties of the eigenfunctions 
4>™k- First of all, we know from (7J Lemma 2.7] that the mean value e(-P) of a polynomial 
p ( x ) = Y!i=m c lPl( x ) can be written as e(P) = c H J™c, where c = (c m , c m+ i, . . . , c n ) T . 
A similar characterization can be found for the variance var(P). 

Lemma 3.1. For a normalized polynomial P(x) = ^IL m QPK 3 -)? we have the following 
characterization of the variance var(P): 

var(P) =c H [J n ] 2 c + b 2 n \c n \ 2 -(c H 3 n c) 2 , if P G U n , 

var(P) = c H [J™] 2 c + b 2 Jc m \ 2 + b 2 n+1 \c n \ 2 - (c H J™c) 2 , if P G II™, m > 1, 

with the coefficient vectors c = (c m , . . . , c n ) T . 



Proof. For m > 1, we denote by p™(x) the vector (p m (x), • • • ,p n (x)) . Then, using the 
three-term recurrence formula (I7J) and the orthonormality relation of the polynomials p\ , 
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3 The localization of the eigenfunctions of P™M X P™ 



V>24,25 (x). 



V&25 ( X ) 



^24,15 (a;)- 



V&15 ( x )- 




Figure 1: Some eigenfunctions ip™ k of the operator P™M X P™ for the Chebyshev poly- 
nomials of first kind. 
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we get for P{x) = E?= m c lPl( x ) G n n , \\P\\w = 1: 

/•l , " 2 

var(P) = / | qa:p; (x) w{x)dx — e(f) 2 

= / I +aipi(x) + b l pi_ 1 (x)) w(x)dx - e(/) 2 

= £ c^p^x) • p n m (x) H J™c w{x)dx + 6^|c m | 2 + b 2 n+1 \c n \ 2 - eU? 

■ H 3™ ( C Pi (x)p~(x)w(x)dx) J™c + b 2 m \c m \ 2 + b 2 n+1 \c n \ 2 - (c H J™c) 



— c [JrT] c + ^ml c m| +& ra +l|c„| — (c J™c) . 

For m = 0, the statement follows analogously but without the term 6^1 c m | 



□ 



Now, we get the following formulas for the expectation value and the variance of the 
eigenfunctions tp. 



m 

n,k' 



Lemma 3.2. For the normalized eigenfunction 1 < k < n — m + 1, corresponding 

to the eigenvalue x™ k , we have 



/ i \ ( i \ 7 2 Pn(x n ,k) 2 

£{Vn,k) = x n , k , var(V'„,fc) = b n+l =^ 



= <fe, var(^ fc ) 



bn+lPn-m{x™k' m ) 2 + 



(19) 
(20) 



Proof. The statements for the mean value e(VCfc) f°U° w directly from the definition of 
the ip™ k as eigenfunctions of the operator P™M X P™ . 

For the variance var(?/>I\) of the normalized eigenfunction ip^u, m > 1, corresponding to 



the eigenvalue x™ k and with the coefficient vector c k given in (17), we can derive from 
Lemma 13.11 that 

var«J = cf [JC] 2 c fe + b 2 m \c m , k \ 2 + b 2 n+1 \c n , k \ 2 - (c£ J™c fc ) 2 

= ( x ™fc) 2 ( c fc^ c A;) 2 + ^ml c m,fc| 2 + I c n,/c P _ ( x ™k c k c k) 2 

62 I 1 2 I t2 I 1 2 



Inserting the coefficients from (17), we get the above result. The same argumentation 
holds also for m = 0. □ 



Remark 3.3. For the case m = 0, the formula (19) for the variance of ip n>k is a special 
case of a variance formula of the Christoffel-Darboux kernel considered in the proof of 
P Theorem 2.2]. 
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If we want the eigenfunction to be localized at the expectation value x™ k , the vari- 
ance of ip™ k should be small, especially if n — m gets large. The question whether the 



variance in ( 19 ) gets small when n is large is linked to a condition known as subexponen- 
tial growth (see [I], |30|). In particular, if the orthonormalization measure w{x)dx is an 
element of the Nevai class M(0, 1), i.e. if the coefficients of the recurrence formula (J6| 
attain the limits lim n _ s>00 a n = and limn^^ b n = |, it is proven in |30| that var(^> nj fc) 
tends to zero as n — > oo. If we restrict the measure w(x)dx to a particular subclass of 
M(0, 1), we can also show in the more general case m > that the variances in Lemma 
13.21 tend to zero as n — > oo. 

Definition 3.4. By M*(0, 1), we denote the set of all measures fx with the following 
properties: 

1. n is in the Nevai class M(0, 1), i.e. lim ra _ 5 . 00 a n = and limn-^oo b n = i, 

2. supp n = [—1,1], 

3 - ^2^=0 \ a n\ + \b n - \\ < OO, 

where a n and b n are the coefficients of the three-term recurrence relation ([7| correspond- 
ing to the measure [i. 



Examples of weight functions lying in the Nevai subclass M*(0, 1) are, for instance, the 
Jacobi weight functions (see |29[ p. 79-81]. 

For a measure u. and the corresponding family of orthonormal polynomials {pi)i^n 1 we 
denote by [i m the orthonormalizing measure of the associated polynomials pi(x,m). In 
particular, the measure /i m is normalized such that fi m ([— 1, 1]) = 1. For a measure fi in 
the Nevai subclass M*(0, 1), we get the following result: 

Lemma 3.5. // u, € M*(0, 1), then also /i m £ M*(0, 1). Moreover, the measures fi m , 
m > 1, are all absolutely continuous on [—1, 1], i.e. d\i m = w m dx. 



Proof. Since the coefficients of the three-term recurrence relation (J7| of the associated 
polynomials pi(x, m) are defined by shifting the corresponding coefficients of the polyno- 
mials pi, the conditions (1) and (3) of Definition 3.4 are obviously satisfied by the measure 
fj, m . The true interval of orthogonality of the sequence of associated polynomials pi(x, m) 
is included in the true interval of orthogonality of the original polynomials pi(x) (see m 
Corollary on page 87]). Therefore, supp fi m C supp/i = [—1, 1]. Since fi m £ M(0, 1) is in 
the Nevai class, also [—1, 1] C supp [i m holds (cf. |29| Chapter 3.3, Lemma 6]) and, thus, 
also the property (2) is satisfied. 

To prove the absolute continuity of fi m we use a result of Nevai |29| Chapter 7, Theorem 
40]. This result implies that if fx G M*(0, 1), then the measure \i consists of an absolutely 
continuous part w{x)dx on [—1, 1] and a point mass a5-\+b8\ on the boundary of [—1, 1]. 
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Hence, it remains to show that for the associated measures fj, m , m > 1 the discrete part 
vanishes. It is enough if we give the proof for the left hand boundary x = — 1. In this 
case, o = is equivalent to the divergence of the sum X]|=oPl( _i i m ) 2 ( c ^- l-^l Theorem 
2.1]). By a technique involving chain sequences, Chihara [3l Formula (2.18)] proved that 
there is a constant C m such that 

K(-l,m + l)| 2 > C m \ Pn+1 (-l,m)\ 2 . (21) 

Hence, by a standard induction argument it follows that X^z=oP'( — ^> m ) 2 ; m — 1 di- 
verges, if ^iLqPi{ — 1, l) 2 diverges. So, to complete the proof we have to show the 
divergence of Y^=oPl{~^-> I) 2 - If A 4 i s continuous at x = —1, then ^2^qPi{ — l) 2 diverges, 
and by ( |21[ ) also YlhLoPl(~ 1 , 1 ) 2 diverges. If has a point mass at x = —1, then by 
another result of Chihara (2J Theorem 3], the measure [i\ cannot have a point mass at 
x = —1. Hence, in this case the sum YlhLoPl(~ 1> I) 2 a l so diverges. □ 

Theorem 3.6. If the weight function w is in the class M*(0, 1), then 

lim var(^ n = 0, lim var(f/>![\,) =0, m £ N, 
uniformly for all k. 



Proof By Lemma 3.5 the measures w m {x)dx lie in the subclass M*(0, 1), hence also in 
the Nevai class M(0, 1). Therefore, by a result of Nevai, Totik and Zhang [30, Theorem 
2.1] we have 

\p n (x,m)\ 2 



lim sup 

n ^°°^e[-i,i] 2^=o b/KH 



0. 



Further, by Lemma 3.5 the associated measures dfi m (x) = w m (x)dx, m > 1, are abso- 
lutely continuous on [—1, 1]. Hence, by [TSJ II, Theorem 2.1], also 

1 



lim 



00 E"=o \Pi{x,m)\ 2 







uniformly on [—1,1]. Therefore, the results of Lemma 3.2 imply that the variances 
v&T(ip nt k) and var(V'™ fc ) converge to zero (independently of the choice of k) as n tends to 
infinity. □ 

Example 3.7. For some particular weight functions w, it is possible to determine the rate 
of convergence of the variance var(ijj n ^) in Theorem 3.6 For instance, if the weight w is 
a generalized Jacobi weight, i.e. if supply - 



w\x) 



-1,1] and 
■1 = ii < t 2 < ■ ■ ■ < U-l < U 



H > -1, 



8=1 
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4 Approximation of localized functions 



then the rate of convergence can be determined as (see |29], Theorem 9.31 and Theorem 
6.3.28) 



var(^ fc ) = fr n+1 n 2 , l<fc<n + l. 

2^,1=0 Pl\ x n,k) n 

So, for generalized Jacobi weights, the convergence of lim™-^ var^^fc) towards zero is 
at least linear. The convergence rate is even faster, if we choose k such that x n ^ is among 
the ./V (iV £ N fixed) smallest or largest roots of p n +\{x). 



4 Approximation of localized functions 



In this paragraph, we are going to investigate how the decomposition of a bandlimited 
function / 6 II™ in the eigenfunctions ij)™u can be used to approximate functions that 
are well-localized at a point or a subinterval of [—1,1]. In this case, not all of the 
eigenfunctions are needed for a good approximation of the function /. We will 

show that mainly only those eigenfunctions are needed that are located themselves in 
the region in which / is concentrated. 

From now on we assume that the weight function w lies in the Nevai subclass M*(0, 1). 
Then, for the Hilbertspace 

L 2 ([-l, l],w)Q II m _i := span{pi : I > m} 

we can introduce an isometric isomorphism S m by 

S m : L 2 ([—l, 1], w) II m _i ->■ L 2 ([-l, l],w m ) : S m pi(x) := pi„ m (x,m), I > m. 

(22) 

If the functions (j) n - mt k, 1 < k < n — m + 1 denote the eigenfunctions of the operator 
P n - m M x P n _ m on the Hilbert space L 2 ([— 1, l],w m ), we can deduce from (18) that 



n—m,k 



» 



holds. Further, for e m > we say that a continuous function f £ L 2 ([—l,l],w) Q II m _i 
is e m -concentrated on an interval A C [—1,1] if 



-l,l]\A 



\S m f (x)\ 2 W m {x)dx < € 2 



2 

w • 



An e m -concentrated function / can be approximated as follows: 

Theorem 4.1. Let f € L 2 {[— 1, 1], w) 0ll m _i be continuous and e m - concentrated on the 
subinterval A C [—1, 1]. Then, 



lim 



k:x m , eA 

n,k 



< e r 



(23) 
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4 Approximation of localized functions 



If A = [cos a, cos (3], the number of eigenvalues x™ k in A is asymptotically given as 

#{k:x™ k £A} a -f3 



lim 

n-s>oo [n — m) 



7T 



Proof. We use the isomorphism S m to shift the error term from the Hilbert space 
L 2 ([-l,l],u>)en m _! onto L 2 ([-l,l],w m ): 



f- E 



£>mf ^ ^ {Srnf)4'n—m,k)w m 4 > n—m,k 



k:x m ,£A 

n,k 



(24) 



For an arbitrary JV 6 N, we can assume without restriction that n is large enough 
such that N < n — m. By P/v = ^2 k=0 (S m f , pk(- , m^wp^^ , m) , we denote the best 
approximation of S m f in the subspace Ily of L 2 ([— 1, l],w rn ), and by 



EN(S m f, w r 



inf ||5 m /-P| 



\S-mf P 



N \\Wn 



the corresponding error term. Now, using (24) and the triangle inequality twice, we get 



/- E </><*>«<» 



P 



A? 



Af> < rn—m.k)w m H , n—m,k 



k:x m ,eA 

n,k 



(25) 



< 



S„J-Pv+ E ( P N-Smf 



i (Pn—m^k) ■w m < Pn—m,k 



< 2E N (S m f, w m ). 



From the spectral Theorem 2.1 we know that the eigenfunctions (j) n - m k are ) U P to a 



normalizing factor, the fundamental polynomials of Lagrange interpolation with respect 
to the nodes x™ k , 1 < k < n — m + 1. In particular, since P/v 6 II jv C II n _ m , we have 
(cf. gB Section 3.4]) 

(PN,4>n-m,k)w m = ^n,k P N(Xn,k)- 

Hence, if we define the bounded function g on [—1,1] by 



P N (x) if xe [-l,l]\A, 
if x e A, 



then the sum 



E ^A r ( x ™fc) K ™A: < An-m,fc 



corresponds precisely to the Lagrange interpolant of 5 at the nodes 1 < k < n—m+1. 
Therefore, by the Erdos-Turan-Theorem (the original result can be found in [9], in our 



14 



4 Approximation of localized functions 



case we need [15 , Chapter 3, Theorem 2.5] with the parameters A n = B n = 0) we get in 
the limit n — > oo: 



lim 

n— >oo 



Pn- e ft 



N ■> < Pn—rn.k/w rn < Pn—m,k 



k:x m ,eA 



lim 

n— >oo 



E ( P N,4>n-i 



m,k/w m < rn—m,k 



g{x) 2 w m {x)dx = I P N (x) 2 w m (x)dx. 
-1 J[-1,1]\A 

Also by the triangle inequality the following estimate holds: 



(26) 



■1,1]\A 



P N (x) 2 w m (x)dx 



-1,1]\A 



S m f(x) 2 w m (x)dx 



< E N (S m f,w r , 



(27) 



Combining (|25|), (|26|) and (|27|), we can conclude for n — > oo: 



lim 

n— ¥co 



f- E </><*>«c 



71, fc 



HUM 



< lim 

n— >oo 



+ lim 

n— >oo 



+ lim 

n— »oo 



/- E 



^- E ^ 



iVi ( Pn—m,k/w m 'Pn—rn,k 



n,fc 



AT; < rn—m,klw rn ( Pn—m,k 



[-i,iM 
< 3E N (S m f,w m ). 

Since iV was choosen arbitrarily, we finally get 



PAr(x) 2 w m (x)(ix 



lim 

n— »oo 



/- E 



S m f(x) 2 w m {x)dx 



0. 



Inequality (23) now follows from the fact that / is 6 m -concentrated on .A. 



Since the weight function uu is in the class M*(0, 1), Lemma 3.5 ensures that also the 
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4 Approximation of localized functions 



associated weight functions w m are in M*(0, 1). This implies suppu> m = [—1,1] and, by 
|29[ Theorem 7.29], that the restricted support of w m on [—1, 1] has measure 2. Therefore, 
by a well-known result of Erdds and Turan (see [8], |T^) w m {x)dx is an arc-sine measure 



which implies the second statement of Theorem 4.1 



□ 



Remark 4.2. The second statement in Theorem 4.1 is not a new result and intended here 
only as an additional information on the asymptotic number of eigenfunctions involved 
in the approximation process. It is a special case of a general property that for a large 
class of orthogonal polynomials the asymptotic distribution of the zeros is given by the 
arc-sine measure. For weights as the functions w m this was proven by Erdos and Turan 
in |10| . Far more general conditions leading to the arc-sine property are elaborated in 
[8]. In particular, it can be shown that every measure in the Nevai class M(0, 1) has this 
property (see \29\ Theorem 5.3]). 



If a polynomial P £ U™ is localized at the end points x = — 1 or x = 1, or if P has a 
small variance var(P), we obtain the following error estimates: 

Theorem 4.3. Let a > and I- and I + denote the Intervals I_ = [—1, —1 + a] and 
I + = [1 — a, 1]. // P G n™, \\P\\w = 1; is localized at the boundary points of [—1, 1], we 
have the following error bounds: 



p- E ^<k) w ^ 



in 
n.k 



n,k 



n, k 



< 



< 



l + e(P) 



1 - e(P) 



a 



(28) 
(29) 



Further, if I = [e(-P) — a,e(P) + a] C [—1, 1] ; we get the following error estimate: 

2 



P~ E ( P ^n,M 



m 
n.k 



n,k 



< 



var(P) 



(30) 



Proof. For P £ 11^, we have 



P 



E 



,m \ 1 2 
J n,k/w\ ■ 



< 



(P,^k) w <k = E K p '^. 

fc*£*e[-i.i]V- 

1 n— m+1 

p,< fc > w | a (i + < fc )<- E K^< fe )-i 2 ( 1 + <* 



fc=l 
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5 An uncertainty principle for the operators M x and P™ 



Since \\P\\l = EkZT +1 m^ 



"" 12 " 1 Ek=T +1 XnA(P^n,kU 2 = <P), we get 



the stated bound (28). In a similar fashion, the bound (29) can be proven. To prove 



(30), we proceed also in a simalar way. 



p- E (P'flbMZ 



k:x m ,£l 



E K p >c 

* : *£*e[-i.i]V 



m \ |2 
fc/u>l 



n—m+1 



fc=l 



<^2 E \(P^nM(e(P)-X™ k ) 2 <- 2 £ |(P,< fc ) M | 2 ( £ (P)-< fe ) 5 

k:x™ k e[-l,l)\I 

n-m+l /px 

= - 2 E i(^^i 2 ((< fc ) 2 -^) 2 )<^. 



fc=l 



□ 



Remark 4.4. Given a normalized polynomial P £ II™, we consider the discrete density 
function p by 



p(x) 



({P^™kU 




if x = x™ k , k = 1, . . . n — m + 1, 
otherwise. 



Then, we can interpret the results of Theorem 4.3 as versions of the Markov and the 
Chebyshev inequality for a />distributed random variable, (cf. |31| p. 114]). 

5 An uncertainty principle for the operators M x and P™ 

We are now going to discuss an uncertainty principle related to the operators M x and 
-P™. In particular, we will discuss the trade off between the space localization of / at 
the boundary points x = 1 and x = —1 of [—1,1] and the frequency localization of / in 
the polynomial subspace II™. The obtained results are very similar to the uncertainty 
principle stated in the theory of Landau, Pollak and Slepian (see |14| . |22j). However, 
the fact that M x is not a projection operator will lead to coarser statements and in 
some extent to differences in the proofs compared to the original setting. A detailed 
proof of the uncertainty principle in the Landau-Pollak-Slepian theory can be found in 
Chapter 2.9] and |22| . An abstract version of the Landau-Pollak-Slepian uncertainty 
principle involving two arbitrary projection operators on a Hilbert space can be found in 
|20[ Part 1, Chapter 3]. An extension of the Landau-Pollak-Slepian uncertainty to more 
general weight functions is given in |26j . 



The main results of this section are summarized in Theorem 5.5 and illustrated in Figure 
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5 An uncertainty principle for the operators M x and P™ 
[2] The proof of the statements in Theorem 5.5 is splitted into four lemmas. We define 



<!■■= iiir/ii« = £i</,p*>«i a 

k=m 

and start with the first auxiliary result. 

Lemma 5.1. Lef f, \\f\\ w = 1, be a fixed normalized function. Then, for every < 
P < C(/) there exists a normalized function g, \\g\\ w = 1, such that e(g) = e(f) and 

CO/) =P- 

Proof We choose k > I > n + 1 big enough such that the three largest eigenvalues xi,x% 
and X3 of the Jacobi matrix 3 l k are larger than e(/). This is possible since the weight 



function w lies in the class M*(0, 1) and Lemma 3.5 ensures that also the associated 
measure wi(x)dx G M*(0, 1) is absolutely continuous on [—1, 1]. Let ipi,ip2 and ^3 denote 
the corresponding eigenfunctions in U k . Further, we define V as the 3-dimensional vector 
space spanned by ipi, ip2, and 1P3, and Py as the orthogonal projection operator from 
L ([— 1, onto V. Now, we take ip as a normalized vector in F that is orthogonal to 

the plane spanned by the vectors Pyf and PyM x f. Then, e(tp) > e(f) and {M x f, ip) w = 
0) (/> V')w) = 0- I n the same way, we construct a normalized vector 99 G IT^ with £(92) < 
e(/) and (xf,(p) w = {f,(f) w = 0. Now, since e(/) is a continuous functional, by the 
intermediate value theorem we can find a normalized polynomial (j> G Ilf with = a 
and {M x f,(J)) w = (f,(f>) w = 0. Then, we define 

S(x) := Vl - A/0) + VX<f>(x), AG [0,1]. 

In this way we get a normalized function g with \\g\\ w = 1, CO?) = (1 ~~ A)7r™(/) and 

e( & ) = l-Ae(/) + Ae(0) = e(/). 

□ 



By x™ min and x™ max , we denote the smallest and the largest root of the associated 
polynomial p n _ m+ i(x,m). Then, we have as a second auxiliary result: 

Lemma 5.2. If 'Cmiii < e(f) < x™ max , then n™(f) can attain all values in the interval 

[0,1]. ' " 

Proof We denote by ^™ max and V'JTmin ^ ne normalized eigenfunctions corresponding to 
the eigenvalues £™ max and x™ min , respectively. Now, for x™ min < a < x™ max , we define 
the function / by 

/ Q _ x m. \ 2 / m _ \ 3 
J n,min | , m I J 'n,max " j < m 



m I ^n,max 1 I m _ T m I rn.min- 
x n,min / \ n,max ^n.min / 
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5 An uncertainty principle for the operators M x and P™ 

Then, tC(/) = H/IU = 1 and 

r(f\= *n,tmn w ^n,max " m 

W / m _ m J 'n,max t m _ m "°n,mm 

x n,max n,min J 'n,max ^jijinm 

Now, Lemma |5 . 1 1 implies the statement. □ 



Lemma 5.3. If £™ max < e(/) < 1, i/ien 7T™(/) can attain all values in the range 
< tt™(/) < • If —1 < < ^JTmin; ^ en n n(f) can attain all values in the 

rar^e < <(/) < j^-. 

Proof We will prove the statement only for the interval [£™ max ;l)> the statement for 
(— l,£™ max ] follows by an analagous argumentation. Since w(x)dx G M*(0,1), we can 

k > / > n + 1 large enough such that 1 — x l k max < e for an 



5.1 



choose as in Lemma 

arbitrary e > 0. Then, for the eigenfunction i/), max G LT^ we have ^^(V^rnax) = an< ^ 

1 > e (Vi, max ) = 4,max > 1 ~ e ' NoW > We define 

5 ( X ) = V^™ ma >) + v 7 !^!,^), A G [0, 1]. 

Then, 

1 — A(l — x™ max ) > e(g) = Ax™ max + (1 — A)x^. max > 1 — e — A(l — x\ 

^ 1 — e — ^(l — ^njinax); 

and TT™(g) = A. Therefore, we get for ir™(g): 

> TTn (5) > 



1 T m 71 w ' 1 T m 

x ■''n.max x ■''r^max 



Since e > can be choosen arbitrarily small, we get the desired result from Lemma 5.1 
□ 

Up to now, we showed that most points (s(f), 7T™(/)) m the rectangle ( — 1,1) x [0,1] 
can be attained for / G L 2 ([— 1, 1], w). However, the next Lemma 5.4 demonstrates that 



tuples (e(f), tt™(/)) in the upper left and right corner of (—1, 1) x [0, 1] are not allowed. 
Lemma 5.4. If x™ max < e(f) < 1, the values of ir™{f) are restricted by 

(e(f) + l)i(x- max + 1)1 + var(/)l(var(/) + (1 + e(f))(s(f) - x™ max ))§ 



<w <- 



var(/) + (e(f) + 1)2 

(31) 
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5 An uncertainty principle for the operators M x and P™ 
For — 1 < e{f) < x™ min , the values ofir™(f) are bounded by 



, (1 - e(f))Hl - < min )2 + var(/)2 (var(/) + (1 - e(f))(e(f) - < min ))^ 
<(f) 2 < — 



var(/) + (l-e(/)) 2 

A simpler but less accurate upper bound for 7r,™(/) is given by 

1 1 



(32) 



< 2 + 2 (e(/)<max + (1 " ^(Z) 2 )^! " Kmax) 2 ) 1 ) (33) 



and 

1 1 



<*(/) < 2 + ^(^/Xmin + (1 - ^(Z) 2 )^! " Kmin) 2 ) 1 ), (34) 

for e(f) in the intervals [x™ max , 1) and (— l,a:^ min ], respectively. 



Proof. We will just prove the inequalities (31) and (33). Inequalities (32) and (|34|l follow 



up to some minor modifications with the same argumentation. Since for 7r™(/) = both 



(31) and (33) are satisfied, we will from now on assume that 7T™(/) > 0. Further, we will 
use the operator Mx+i on L 2 ([— 1, X],w) defined by Mx+i f(x) := -^f(x). 

For a normalized function / £ L 2 ([— 1, 1], w) the two functions g\ = ^ M ^ ^ Mx+i f 

and g2 = \\pmf\\ ^nf are a l so normalized. The sum of the angular distances between 
the vectors g\ and /, and gi and / is always larger than the angular distance between 
gi and g 2 , i.e. 

arccosRe((7i, f) w + arccosRe(g2, f)w > &rccos Re (<?i, (72)10 ■ (35) 

i i 

We define the positive selfadjoint operator M% +1 by MJ +1 := M nn^- Then, for the 



term Re(<7i, <?2)iu> we can find an upper bound using the Cauchy-Schwarz-inequality and 
Theorem O 



Re(5i, 52)10 < 1(51,52)^! 



|M^i/|U||P-/||, 



KM^f^^P^f)^ ^ A /(M 2 ± 1 /,/) WA /{M £ ±iP-/,P-/), 



< 



llM^+i/IUHP^/IU " llMs+i/yiP™/^ 
2 2 

^(M^fJ) w ^ X -^f^y/(P™f,P™f) w ^(M^fJ) w ^/^f^ 
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5 An uncertainty principle for the operators M x and Pj 



Now, if we rewrite the expressions (Mx+i f, /),„ and llMg+i f\\ w in terms of e(f) and 

2 2 

var(/), we get 



yW(/) + Rjf+W ' 
£(/) + 1 



Re(g 1 ,f) w = , 

v / var(/) + (e(/) + l) 2 

Re(g 2 ,f) w = ^™f. 



Inserting this into inequality ([35]) , we obtain 



arccos ■ 



+ arccos y 7r„ / > arccos — . (ob) 



VW(7) + (e(f) + 1)2 v nJ " ' yW(/) + ( e (/) + 1)2 

Applying the cosine addition formula, this inequality can be rewritten as 

>(/) + l)^«max + l) 



var(/) + (£(/) + 1)2 



/ (£(/) + 1) 2 \i/ (£(/) + i)« max + i; 

+ l var(/) + (e(/) + 1)2 J I var ( /) + ( £ (/) + 1)2 

+ l)i« max + 1)5 + var(/)l(var(/) + (1 + e(f))(e(f) - < max ))^ 
var(/) + (s(f) + l) 2 



Hence, inequality (31) is shown. To prove inequality (33), we consider inequality (36). 



For < a < b < 1, the function arccos bt — arccos at is a decreasing function of the 

. I ll n ora f™«, ;f „r Q co + ^ — . / Z ",max + 1 ^ 7 _ / e(/) + l 



variable t G ±]. Therefore, if we set a = y "- m 2 ax <b=yj < 1 and 

/ 2( £ (/) + l) x§ = / (M^iUU n| 
Vvar(/) + (£(/) + 1)2 v 1 ^(M a ^/,Af.±i/)J ~ ' 



we get in inequality (36) the upper bound 



/ e (/) + l\2 / , / ^"max + 1 \ 2 

arccos I I + arccos y 7r™ j > arccos 



or equivalently 



y^\7 < + l)5(x- max + 1)3 + (1 - - < max )s). (37) 



Taking the square of both sides in (37), we obtain precisely inequality (33). □ 
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Now, we introduce the functions 71 (x) and 72(21) by 

1 1 



71 (x) :[< max ,l)^R: 7l (x) := - + - (x< max + (1 - x 2 )? (1 - (x 
72(2?) :(-l,< min ) -> R : 72 («) : = 3 + ^ ( x <min + (1 - ^(1 - (x™ min ) 2 )2). 
and the following subdomains of the rectangle (—1, 1) x [0, 1] (see Figure [2]): 

A := {(x, y) G (-1, 1) x [0, 1] : y < j^g^, y < T+^g-} U {(x™ max , 1), (x™ min , 1)}, 
Bi := {(x, y) G (x™ max , 1) x [0, 1] : y > , y < -yi(x)}, 

ri .max 

B 2 :={(x,y) e (-l,x™ min ) x [0, 1] : y > ]+ x , y < 72 (x)}, 

Ci := {(x, y) G (x™ max , 1) X [0, 1] : y > 71(a)}, 
C 2 := G (-l,x^ min ) X [0, 1] : y > 72 (x)}. 



Finally, we can summarize the results of Lemma |5.1| |5.2| |5.3| and 5.4 as follows. 



Theorem 5.5. For normalized functions f G L 2 {[— 1, 1], w), all points (e(/), 7r™(/)) in 
the domain A can be attained. All points (e(/), 7r™(/)) in the corners C\ and C 2 cannot 
be attained. 



Remark 5.6. Theorem 5.5 and its proof based on the Lemmas formulated before are 



highly inspirated by the uncertainty relation of the original Landau-Pollak-Slepian theory 



as described in [5j Chapter 2.9], |20| Part 1, Chapter 3] and |22j . Lemma 5.1 reproduces 



statement F in |20l Part 1, Section 3.1, p. 95]. However, since M x is not a projection 



operator, the proof is altered considerably. Lemma 5.3 is an adaption of Case 2 in the 
proof of |22[ Theorem 2]. The idea for the proof of Lemma 5.4 is taken from [20, Part 
1, Section 3.1 E), p. 95] and the proof of Case 3 in |22[ Theorem 2]. Due to the fact, 
that M x is not a projection operator also here the proof differs from the original one. 
Moreover, the resulting inequalities can not be shown to be sharp. Bounds from below 



are given in Lemma 5.3 but it is not yet clear to which extent points (e(f), 7rJ^(/)) can 



be attained in the domains B\ and S 2 . 
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